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We determine the capillary attraction and equilibrium configurations of particles trapped at a 
liquid-fluid interface due to the pinning of their contact-line. We calculate analytically the asymp- 
totic interaction energy between two particles and, numerically, the multi-body energy landscape 
for up to four contacting particles. Our results are consistent with recent experiments. We show 
that a system composed of a large number of such particles behaves as a jammed system. 

PACS numbers: 68.03.Cd, 82.70.Dd, 83.80.Hj 



Capillary interactions arise whenever particles that are 
floating between a liquid and a fluid distort the inter- 
face. These interactions originate from the dependence 
of the energy stored in the interface on the distance be- 
tween the particles. More generally, interactions medi- 
ated by the elasticity of the medium in which particles are 
embedded occur in, e.g., phase-separating binary mix- 
tures [[[) , critical fluids p , bilayer membranes g, ^] , and 
liquid crystals ^ [|. Capillary interactions play signif- 
icant roles in several important technologies 0, includ- 
ing ore flotation || , foam stabilization ||] , and solid sur- 
face nanostructuring by deposition of two dimensional 
colloidal crystal or aggregates (l0[ [TlJ [L2| . 

The interaction between floating objects is well under- 
stood only when the contact-lines can freely adjust their 
shape at the particles' surfaces: the weight of the parti- 
cles induces an interfacial distortion that relaxes by cap- 
illarity, thereby mediating the interaction. Exact numer- 
ical results are only available, however, for the simplest 
case of two parallel cylinders Various approxima- 

tions, e.g., Nicolson's superposition approximation p3, 
have been devised for spheres and cylinders jl5|, [l6|, p"7f . 

Recently, strong capillary aggregations of floating col- 
loidal spheres have been reported in a regime where the 
above gravity-induced interaction is negligible, due to the 
smallness of the particles |Q. This effect has been at- 
tributed to a pinning of the contact-line at the parti- 
cles' surface. The corresponding long-range interfacial 
distortion has a quadrupolar shape, which mediates an 
interaction very strong in units of k^,T. Similar be- 
haviors were recently observed in experiments involving 
photolithography-fabricated curved discs having an hy- 
drophobic and an hydrophilic side. It is then the edge of 
the curved disc that effectively pins the contact- line [^9) . 

In this paper, we model both experiments in an uni- 
fied way and we produce exact results for the interaction 
between such "trapped" particles. We obtain the asymp- 



totic interaction analytically without using uncontrolled 
superposition approximations. We then compute numer- 
ically the exact interaction between two particles, and 
the many-body interaction between three and four par- 
ticles, in order to understand collective behaviors. We 
finally attempt to connect our results with the jammed- 
fluidized transition for attractive particles under external 
stress ||(], ^l) , thereby proposing these systems as models 
for micro-rheological studies. 

For definiteness, we consider floating particles that 
posses a fixed, non-planar line that pins the fluid's in- 
terface (as in the above examples When such 
particles interact, their contact-lines rise, tilt and ro- 
tate rigidly with the particles; however, each line re- 
mains fixed in the local frame attached to its particle. 
For simplicity, and in agreement with Ref. we as- 
sume that the contact-lines depart weakly from a cir- 
cle, thus excluding elongated particles. We also assume 
weak distortions of the interface height h(x, y) relative 
to the plane (x,y). This implies that the projection of 
each contact-line onto this plane can be approximated 
by a circle, even when the particle is tilted. We de- 
scribe the i-th particle's contact-line by the height func- 
tion Ci(4>i) ( see Fig-0)j that we expand in Fourier har- 
monics: C<(0») = Co; + H„[C m cos(n^,) + S ni sin(n^i)]. 
Due to the fixed shape of the contact-line in the parti- 
cle's frame, some of the above coefficients are constrained 
while others are free. Co;, which describes how the par- 
ticle rises, is free (Fig.|l|). The dipolar coefficients Cu 
and Su, which parametrize the tilts of the contact-line 
relative to the z-axis, are also free. The lowest-order 
multipole bearing the information concerning the shape 
of the contact-line is therefore the quadrupole We 
assume explicitely 



CiOfc) = C 0l + Cu cos^j + Su sin fa 

+ H 2l cos [2 (<f>i - 0i)] + H 4i cos [4 i 



•&)]■ (i) 
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FIG. 1: Representation of the contact-line of the i-th particle 
and of its projection onto the reference plane. 



The quadrupolar coefficient Hn is a constant, associated 
with the shape of the contact-line, and Pi is the free vari- 
able that represents the rotation of the particle about its 
normal n (Fig.|l|). For the sake of completeness we have 
added the lowest-order multipole that does not break 
mirror-symmetry, if 4$ cos [4(0* — Pi)]. The octupolar co- 
efficient H^i is a constant that describes how the contact- 
lines departs from a purely quadrupolar shape. The two 
points where the height modulation of the contact-line 
relative to its average plane is the largest define an axis, 
evidenced by the dashed line in the figures. 

For small deformations, the excess free-energy associ- 
ated with the interface's surface tension 7 is given by 



F = 



d 2 r {Vhf 



(2) 



where r = (x,y). For micron-sized particles as we con- 
sider here, gravity is negligible [jTsfl . At equilibrium the 
first variation of F vanishes, which yields (see Fig. fil) 



d 2 r (V 2 /i) dh 



0. 



(3) 



where dsi = aidfa. Thus V 2 /i = 0. For fixed projected 
positions and rotations Pi of the particles, reduces 
to SCoi + 8C\i cos0j + 8S\% sin^i, since the /Vs and the 
coefficients Fin and in Eq. (|l|) are constant. Hence 
the boundary equilibrium (freedom of the particles to 
rise and tilt) is expressed by the vanishing of the zeroth 
and first-order Fourier harmonics of dh/dri on T>i, which 
implies that the expansion of dh/dri starts at n — 2: 



dh 

dr. 



{r. l = a l , fc) =} c ni cos (n^i) + s ni sin (n<f>i) . (4) 



?i=2 



The condition V 2 h — implies the integral equation: 



Mr') = W 

i J Vi 



dri 



dri 



dsi, (5) 



known as the second Green identity p^|, where G(r, r') = 
log |r — y'\/2it is V 2 's Green function. Evaluating Eq. (||) 
for r' on the ZVs with the Fourier series (0) and (g), and 
with the multipolar expansion: 



G(r,r') = - 

Z7T 



v--> r n 

logr'-^J— cos [n(cf> - (/>')} 



n=l 



(6) 



we obtain a linear set of equations for the unknown vari- 
ables Co%, Cu, Su, and c ni , s ni for n > 2. Given Pi, H 2l 
and Ha, these variables determine completely both h 
and dh/dri on the Pi's. Finally, Eq. (||) determines h(r) 
everywhere. The corresponding free-energy is obtained 
by integrating Eq. (0) by parts and using V 2 /i = 0: 



F = 1 Jd 2 rV-(hVh) = ~lj2 a >J Q 



dh 

dri 



= -— ^2 ai[H ki (c k iCos{kPi) + Skism(kPi))]. (7) 

i; fc=2,4 

Two particles. — Consider two particles (i — I, r) sep- 
arated by a distance d (Fig. ||). Solving analytically 
the linear system determining the unknown profile co- 
efficients for d -> 00 yields F(d) = F + F 4 (d) + 0(d- 6 ), 
where the leading-order interaction is 

2 2 

F 4 = -127r 7 ^-ff 2£j ff 2r cos [2 (p e + P r )} . (8) 

It coincides with the result found in Ref. pl| from a su- 
perposition approximation in the spirit of Nicolson |L4| . 
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FIG. 3: Normalized interaction energy AF/fH2 2 for three 
identical particles with H4 = 0.25 H2 as a function of the 
configuration angle ijj. The two branches come from the min- 
imization with respect to the angles ft . The barrier associated 
with the minimum A is ~3.6 x 10 4 ksT for H2 = 0.1 /Ltm. 

Note that it does not depend on the H^s. At fixed d, 
the particles interaction is minimized when the particles' 
axes are symmetrically rotated with respect to the line 
joining them (ft = —ft), as if the particles were geared. 
Then, the force acting between them is attractive, asymp- 
totically given by 

, . a,f^a T ^ HopHor , , 

/00(d) = ~48tt 7 ^ 2r . (9) 

The exact interaction energy AF(ci, ft,ft) can be ob- 
tained for any distance by solving numerically the set of 
equations for the profile coefficients after truncating the 
Fourier expansions at a large but finite order, and check- 
ing for convergence, /oo turns to be a very good approx- 
imation for a wide range of separations, as evidenced in 
Fig. ||, for two identical particles with Hu = 0. The 
corresponding force at contact is however ~65% weaker 
than /oo. The associated interaction energy is —2.0 7H2 2 , 
~20% weaker than the approximation given by F4.. For 
7 ~ 70 erg/cm 2 (water), a ~ 1 /im, and FL-ija ~ 0.1, this 
gives ~ -1.4 x 10~ 8 erg (~ -3 x 10 5 k B T). Note that 
the contact force actually depends on the details of the 
higher Fourier components. 

When the H^s are zero (purely quadrupolar parti- 
cles), the numerically determined interaction turns out 
to be proportional to cos [2 (ft + j3 r )\ at any separa- 
tion: the energy minimum is degenerate, corresponding 
to ft = —f3 r - Physically, this "gearing" condition is fa- 
vored because it allows the heights £i of the particles' 
contact-lines to match, thereby preventing a strong local 
interfacial distortion jl9| . Our exact calculations further 
showed that the presence of nonzero octupolar correc- 
tions Ha actually lifts this degeneracy: the minimum 
energy is attained when ft = ft — (particles' axes 
parallel to the separation vector) , as intuited in Ref . p9| . 




FIG. 4: Normalized interaction energy AF/7H2 2 for four 
identical particles with H4 — OAH2 as a function of the an- 
gle ip r specifying the position of the rightmost particle. 

Indeed, it is best to match the heights where the contact- 
lines have a sharper distortion, i.e., at the tips of the axes. 

Three particles. — Owing to the strength of the above 
attraction (~ 10 5 k^T), we study the equilibrium states of 
three particles by assuming: i) that each particle contacts 
at least another particle, ii) that the particles can rotate 
around each other, in agreement with the experiments 
of Ref. |ist| . We consider generic particles, i.e., H 4i ^ 0. 
Without loss of generality, we fix two particles at contact, 
and we explore the energy landscape as a function of the 
position of the third particle, described by the angle i/j 
(Fig. ||). For each tp we minimize the energy with respect 
to the rotations ft of the three particles. 

We obtain two distinct branches exhibiting three 
metastable minima (see Fig. ||). The latter almost per- 
fectly satisfy the above-mentioned "gearing" rule. The 
curved-\ine configuration A suggests a three-body effect, 
since the particle's axes are not parallel to one-anothcr, 
contrary to the pairwise case. Indeed, repeating the cal- 
culation with pairwise interactions instead of the exact 
many-body interaction yields a straight line. The most 
stable, packed configuration is the triangle C with a ra- 
dial orientation of the particles. One might argue that 
it best relaxes the strong interfacial distortion located at 
the facing tips of the particles' axes. The orthoradial 
triangle B is also very stable. 

Four particles. — To determine the equilibrium states 
of four contacting particles, we explore the energy land- 
scape by fixing the positions of three particles through 
the angle t() r , and minimizing with respect to both the 
position tp£ of the fourth one and all the rotation angles 
ft (see Fig. ||). We obtain three branches and six equi- 
librium states. 

As a new indication of a many-body effect, the "gear- 
ing" rule is strongly violated in configuration D and 
slightly violated in C. Contrary to the three-particle case, 
the orthoradial square B is more stable than the radial 
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mal systems proposed in Ref. pl| , with a characteristic 
power-law a y ~ 

It should thus be very interesting to experimentally 
study the micro- rheological behavior of such particles' 
aggregates. Indeed, direct observations would be possible 
because the system is two-dimensional; furthermore, as 
we have shown, this is one of the few systems where exact 
many-body calculations are possible. 

We thank D. Johannsmann, A. Ajdari, E. Raphael, 
D. T. Wu, and P. Flaud for interesting discussions and 
comments. 



FIG. 5: Tentative jamming phase diagram extrapolated from 
the results for four particles. The dashed line is the power-law 

CFy ~ 



square E. One might suspect that since the gap between 
the particles is larger than in the triangular case, it is 
more important to relax the distortion stored at the axes' 
tips on the sides of the square rather than in its cen- 
ter. Note however that such reasonings become more 
and more delicate as the number of particles increases. 
The fact that B and E are the most stable configurations 
suggests, as proposed in Ref. that the ground state 
for many particles should be a square lattice with B, or 
equivalently E, as unit cell. We find again the curved line 
A and different triangular arrangements with a fourth 
attached particle (C, D, F). Our results do successfully 
reproduce the various structures that can be locally seen 
in the experimental pictures of Ref. Jp| . 

Such floating particles can be considered as a model 
system to study jamming phenomena in complex flu- 
ids. Jamming occurs in the presence of either repulsive 
or attractive forces |21j| , when a disordered system gets 
trapped in phase space, for kinetic or energetical reasons, 
and thus exhibits a yield stress . Aggregates of a large 
number of the "trapped" particles discussed in this pa- 
per should form by accretion of individual particles or 
small clusters. In the absence of irreversible sticking, at 
each accretion the particles will rearrange until they fall 
in one of the many metastable configurations. They will 
thereby form a disordered structure with a yield stress. 
Such colloidal aggregates are therefore more akin to a 
foam |2^| than to a collection of sticky spheres j24 ]. 

To test the relevance of this description in terms 
of jamming, we estimate the dependence of the yield 
stress <7y on the volume fraction (/>, from the energetical 
landscape for four particles (Fig. ^). For each metastable 
state, we set (f> oc 1/L 2 , where L is the distance between 
the farthest particles, and we estimate a y ~ E^/aL, 
where E\, is the energy barrier towards the nearest mini- 
mum, and ~ Eb /a the force required to rearrange a given 
state. This yields the diagram of Fig. |[ which is quite 
reminiscent of the phase diagram for attractive ather- 
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